*******************************************************
*******************************************************
*******************************************************
*** CREATED BY JONATHAN KING AND JESSICA SCHOENHERR ***
***** REPLICATION DATA FOR "A MATTER OF OPINION?" *****
*******************************************************
*******************************************************
*******************************************************


use "DeathPenaltyData20220224.dta"

**********************
*** BASELINE MODEL ***
**********************

reg decisionTherm ib10.treatmentGroup
est store dpBaseline

est restore dpBaseline
margins, at(treatmentGroup=(6(1)10))

est restore dpBaseline
margins, at(treatmentGroup=(10)) ///
	at(treatmentGroup=(6)) post
test _b[1._at] = _b[2._at]
* not sig

est restore dpBaseline
margins, at(treatmentGroup=(10)) ///
	at(treatmentGroup=(7)) post
test _b[1._at] = _b[2._at]
* not sig

est restore dpBaseline
margins, at(treatmentGroup=(10)) ///
	at(treatmentGroup=(8)) post
test _b[1._at] = _b[2._at]

est restore dpBaseline
margins, at(treatmentGroup=(10)) ///
	at(treatmentGroup=(9)) post
test _b[1._at] = _b[2._at]

est restore dpBaseline
margins, at(treatmentGroup=(6)) ///
	at(treatmentGroup=(7)) post
test _b[1._at] = _b[2._at]

est restore dpBaseline
margins, at(treatmentGroup=(6)) ///
	at(treatmentGroup=(8)) post
test _b[1._at] = _b[2._at]

est restore dpBaseline
margins, at(treatmentGroup=(6)) ///
	at(treatmentGroup=(9)) post
test _b[1._at] = _b[2._at]

est restore dpBaseline
margins, at(treatmentGroup=(7)) ///
	at(treatmentGroup=(8)) post
test _b[1._at] = _b[2._at]

est restore dpBaseline
margins, at(treatmentGroup=(7)) ///
	at(treatmentGroup=(9)) post
test _b[1._at] = _b[2._at]

est restore dpBaseline
margins, at(treatmentGroup=(8)) ///
	at(treatmentGroup=(9)) post
test _b[1._at] = _b[2._at]

********************
*** PARTISANSHIP ***
********************

* dem 1
* ind 2
* rep 3

reg decisionTherm ib10.treatmentGroup##ib1.partisanship
est store dpPartisanship

est restore dpPartisanship
margins, at(treatmentGroup=(6(1)10) partisanship=(1(1)3))

* democrats across treatment
*compared to group 6
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(7) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 7
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(7) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(1)) ///
	at(treatmentGroup=(7) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(1)) ///
	at(treatmentGroup=(7) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 8
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(8) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(1)) ///
	at(treatmentGroup=(8) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 9
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(9) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

* republicans across treatment
* compared to group 6
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(3)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(3)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(3)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(7) partisanship=(3)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 7
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(3)) ///
	at(treatmentGroup=(7) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(3)) ///
	at(treatmentGroup=(7) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(3)) ///
	at(treatmentGroup=(7) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

*compared to group 8
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(3)) ///
	at(treatmentGroup=(8) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not, p = 0.0529

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(3)) ///
	at(treatmentGroup=(8) partisanship=(3)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 9
est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(3)) ///
	at(treatmentGroup=(9) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

* dems vs. reps within a treatment
est restore dpPartisanship
margins, at(treatmentGroup=(6) partisanship=(1)) ///
	at(treatmentGroup=(6) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0151

est restore dpPartisanship
margins, at(treatmentGroup=(7) partisanship=(1)) ///
	at(treatmentGroup=(7) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0299

est restore dpPartisanship
margins, at(treatmentGroup=(8) partisanship=(1)) ///
	at(treatmentGroup=(8) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0260

est restore dpPartisanship
margins, at(treatmentGroup=(9) partisanship=(1)) ///
	at(treatmentGroup=(9) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0001

est restore dpPartisanship
margins, at(treatmentGroup=(10) partisanship=(1)) ///
	at(treatmentGroup=(10) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0001

**************
*** GENDER ***
**************

reg decisionTherm ib10.treatmentGroup##i.female
est store dpGender

est restore dpGender
margins, at(treatmentGroup=(6(1)10) female=(0(1)1))

* men across treatment
* compared to group 1
est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(6) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(0)) ///
	at(treatmentGroup=(6) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(8) female=(0)) ///
	at(treatmentGroup=(6) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(7) female=(0)) ///
	at(treatmentGroup=(6) female=(0)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 7
est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(7) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(0)) ///
	at(treatmentGroup=(7) female=(0)) post
test _b[1._at] = _b[2._at]
* p < .10; 71.19 (7) vs. 62.94 (6)

est restore dpGender
margins, at(treatmentGroup=(8) female=(0)) ///
	at(treatmentGroup=(7) female=(0)) post
test _b[1._at] = _b[2._at]
* yes; p=0.0072; 71.19 (7) vs. 57.82 (7)

*compared to group 8
est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(8) female=(0)) post
test _b[1._at] = _b[2._at]
* no, p = 0.0640

est restore dpGender
margins, at(treatmentGroup=(9) female=(0)) ///
	at(treatmentGroup=(8) female=(0)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 9
est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(9) female=(0)) post
test _b[1._at] = _b[2._at]
* no

* women across treatment
*compared to group 6
est restore dpGender
margins, at(treatmentGroup=(10) female=(1)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(1)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(8) female=(1)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(7) female=(1)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* yes (p<0.1); 60.78 (6) vs. 53.66 (7)

*compared to group 7
est restore dpGender
margins, at(treatmentGroup=(10) female=(1)) ///
	at(treatmentGroup=(7) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(1)) ///
	at(treatmentGroup=(7) female=(1)) post
test _b[1._at] = _b[2._at]
*no

est restore dpGender
margins, at(treatmentGroup=(8) female=(1)) ///
	at(treatmentGroup=(7) female=(1)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 8
est restore dpGender
margins, at(treatmentGroup=(10) female=(1)) ///
	at(treatmentGroup=(8) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(1)) ///
	at(treatmentGroup=(8) female=(1)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 9
est restore dpGender
margins, at(treatmentGroup=(10) female=(1)) ///
	at(treatmentGroup=(9) female=(1)) post
test _b[1._at] = _b[2._at]
* no

* men vs women, within a treatment
est restore dpGender
margins, at(treatmentGroup=(6) female=(0)) ///
	at(treatmentGroup=(6) female=(1)) post
test _b[1._at] = _b[2._at]
* no 

est restore dpGender
margins, at(treatmentGroup=(7) female=(0)) ///
	at(treatmentGroup=(7) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0002

est restore dpGender
margins, at(treatmentGroup=(8) female=(0)) ///
	at(treatmentGroup=(8) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpGender
margins, at(treatmentGroup=(9) female=(0)) ///
	at(treatmentGroup=(9) female=(1)) post
test _b[1._at] = _b[2._at]
* no, p = 0.0859

est restore dpGender
margins, at(treatmentGroup=(10) female=(0)) ///
	at(treatmentGroup=(10) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 0.0858

******************
*** FULL MODEL ***
******************

reg decisionTherm ib10.treatmentGroup##ib1.partisanship##i.female
est store dpFull

est restore dpFull
margins, at(treatmentGroup=(6(1)10) partisanship=(1(1)3) female=(0(1)1))

* democrat men across treatment
*compared to group 6
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0285 (increase)

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*yes; p<0.05; 69.21 (6) vs. 52.6 (8)

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0262 (increase)

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*yes; p=0.0234; 68.93 (7) vs 52.6 (8)

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no

* democrat women across treatment
*compared to group 6
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no 

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

* Republican men across treatment
*compared to group 6
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0364 (decrease)

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes; p=0.0431; 62.96 (6) vs. 81.2 (9)

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0160

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*yes; p=0.0187; 64.91 (8) vs. 81.2 (9)

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

* Republican women across treatment
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes; p=0.0780; 74.48 (6) vs. 62.68 (9)

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*yes; p=0.0494; 74.48 (6) vs. 60.85 (8)

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*yes; p=0.0296; 74.48 (6) vs. 59.08 (7)

*compared to group 7
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*no

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 8
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*no

*compared to group 9
est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

* democrat men vs women, within a treatment
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0190

est restore dpFull
margins, dydx(female) at(treatmentGroup=(6) partisanship=(1))

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0010

est restore dpFull
margins, dydx(female) at(treatmentGroup=(7) partisanship=(1))

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(8) partisanship=(1))

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(9) partisanship=(1))

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(10) partisanship=(1))

* republican men vs women, within a treatment
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(6) partisanship=(3))

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0418

est restore dpFull
margins, dydx(female) at(treatmentGroup=(7) partisanship=(3))

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(female) at(treatmentGroup=(8) partisanship=(3))

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0080

est restore dpFull
margins, dydx(female) at(treatmentGroup=(9) partisanship=(3))

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0260

est restore dpFull
margins, dydx(female) at(treatmentGroup=(10) partisanship=(3))

* men, republican vs. democrat, within a treatment
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(6) female=(0))

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(7) female=(0))

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no, p = 0.0754

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(8) female=(0))

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0010

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(9) female=(0))

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(10) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0001

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(10) female=(0))

* women, republican vs. democrat, within a treatment
est restore dpFull
margins, at(treatmentGroup=(6) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(6) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0008

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(6) female=(1))

est restore dpFull
margins, at(treatmentGroup=(7) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(7) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no, p = 0.0911

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(7) female=(1))

est restore dpFull
margins, at(treatmentGroup=(8) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(8) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(8) female=(1))

est restore dpFull
margins, at(treatmentGroup=(9) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(9) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, 0.0299

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(9) female=(1))

est restore dpFull
margins, at(treatmentGroup=(10) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(10) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, 0.0381

est restore dpFull
margins, dydx(partisanship) at(treatmentGroup=(10) female=(1))



